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ABSTRACT 

In  this  article  we  present  a  multiple  scattering  analysis  of  the  coherent  wave 
propagation  through  an  inhomegeneous  medium  consisting  of  either  random  or 
periodic  distribution  of  scatterers  of  arbitrary  shape.  Both  specific  and 
random  orientations  of  the  scatterers  are  considered.  The  mathematical  unity 
inherently  present  in  the  T-matrix  formalism  for  the  three  wave  fields,  namely 
acoustic,  electromagnetic  and  elastic,  is  employed  in  conjunction  with  suitable 
averaging  procedures  to  formulate  a  self-consistent  multiple  scattering  theory. 
For  a  random  distribution  of  scatterers  we  use  a  configurational  averaging 
procedure,  while  for  a  periodic  distribution,  we  use  a  suitable  lattice  sum 
based  on  crystallographic  theory.  The  information  about  the  orientation  of  the 
scatterers  has  been  incorporated  into  the  T-matrix  of  the  scatterer  itself  thus 
making  formalism  a  convenient  computational  scheme  to  study  the  anisotropic 
effects  in  an  inhomegeneous  medium.  The  statistically  averaged  equations 
obtained  by  the  analysis  are  then  solved  by  using  Lax's  quasicrystalline  appro¬ 
ximation  to  obtain  the  bulk  or  effective  properties  of  the  medium.  Numerical 
results  are  presented  for  propagation  speeds,  attenuation  and  frequency  dependent 
elastic  properties  for  a  range  of  frequencies  to  demonstrate  the  broad  appli¬ 
cability  of  the  T-matrix  method. 

I.  INTRODUCTION 


The  subject  of  wave  propagation  and  scattering  in  inhomogeneous  media  has  become 
increasingly  important  in  many  fields  of  engineering  and  science,  for  e.g.,  in 
studies  of  bubbles  in  a  fluid,  distribution  of  flaws  in  solids  (NDE) ,  ionospheric 
irregularities,  geophysical  exploration,  artificial  dielectrics,  millimeter  wave 
propagation  in  ocean  fogs  and  mists,  and  through  rain,  cloud  and  hail  particles, 
underwater  signal  transmission,  porous  media,  composites,  attenuation  of  waves 
in  biological  media,  remote  sensing,  etc.  All  of  these  problems  are  characterized 
by  a  suitable  statistical  description  of  the  media.  The  similarity  in  statistical 
descriptions  and  the  mathematical  unity  present  in  the  description  of  all  three 
classical  wave  fielis,  namely  acoustic,  electromagnetic  and  elastic,  may  be  taken 
into  consideration  in  formulating  a  unifying  approach  to  all  these  problems.  In 
this  article,  we  present  one  such  unifying  approach  based  on  the  T-matrix  or  the 
null  field  method  that  has  been  recently  formulated  for  all  three  fields  and 
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promises  to  be  an  efficient  computational  scheme  to  analyze  the  scattering  of 
waves  from  several  different  geometries. 

At  any  point  in  a  random  medium,  the  total  wave  fields  can  be  considered  as  a 
sum  of  two  components,  viz,  a  coherent  or  average  wave  (which  is  the  statistical 
average  over  all  possible  configurations  of  the  scatterers  with  regard  to  location 
and  state)  and  an  incoherent  component  due  to  changes  in  scatterer  positions  and 
states  from  configuration  to  configuration.  The  averages  of  the  square  of  magnitude 
of  the  coherent  and  incoherent  fields  are  called  the  coherent  and  incoherent  inten¬ 
sities,  respectively.  The  total  intensity  is  equal  to  the  sum  of  coherent  intensity 
and  incoherent  intensity.  For  a  plane  wave  incident  on  a  medium  containing  random 
particles,  the  coherent  intensity  attenuates  due  to  scattering  and  absorption. 
Incoherent;  scattering  effects  introduce  'noise'  into  the  system  and  cause  fluc¬ 
tuations  in  the  coherent  received  amplitude  and  phase.  In  radar  meteorological 
applications  it  is  important  to  assess  the  incoherent  scattered  intensity  relative 
to  the  total  intensity  in  order  to  relate  theoretical  and  experimental  results 
(e.g.,  power  law  relations  between  attenuation  and  rainfall  rate).  Propagation 
of  the  coherent  wave  is  generally  expressed  in  terms  of  a  bulk  propagation 
coefficient  characterizing  the  scatterer  filled  medium.  Incoherent  effects  are 
usually  determined  by  solving  'approximate'  integral  equations  or  by  solving 
special  forms  of  the  radiative  transfer  equations.  Such  formulations  are  generally 
valid  under  conditions  of  sparse  concentrations  and/or  weak  multiple  scattering 
for  either  Rayleigh  scatterers  or  large  scatterers  which  scatter  primarily  in  the 
forward  direction.  Incoherent  scattering  is  beyond  the  scope  of  this  present 
paper,  and  those  who  are  interested  in  such  analysis  may  refer  to  Ishimaru  (1) . 
Here,  our  attention  focuses  on  formulating  a  multiple  scattering  theory  for 
studying  coherent  wave  propagation  through  an  inhomogeneous  medium  consisting  of 
either  random  or  periodic  distribution  of  scatterers  of  arbitrary  shape. 

Several  theories  and  models  on  diffraction  and  scattering  of  waves  have  been 
pursued  ever  since  Rayleigh's  (2)  analytical  treatment  of  scattering  of  waves  by 
randomly  distributed  particle  to  study  the  color  of  the  sky.  We  cite  here  papers 
that  are  related  our  present  analysis  (3-21).  The  limitations,  difficulties 
and  advantages  of  various  approaches  are  discussed  in  a  recent  review  article  (3) . 
We  also  refer  to  (1)  and  the  references  cited  therein.  Since  scattering  theory 
starts  with  discrete  ensemble  of  inhomogeneities  before  statistical  averaging  is 
carried  out,  the  specific  geometry  and  orientation  of  each  scatterer  can  be 
easily  incorporated  into  the  formulation.  This  is  not  possible  with  continuum 
theories.  The  additional  advantage  of  the  scattering  theory  approach  is  that  it 
is  the  counterpart  of  actual  experiments  performed  for  coherent  wave  propagation 
in  heterogeneous  media  using  elastic,  electromagnetic  or  acoustic  waves  as 
appropriate. 

In  multiple  scattering  theories,  approximations  are  usually  made  for  the  geometry 
and  size  of  the  scatterer  relative  to  the  wavelength  of  incident  waves,  and  the 
distribution  of  scatterers  in  the  medium.  The  api roximations  with  respect  to 
geometry  and  size  of  the  scatterers  are  related.  If  the  scatterers  are  small 
compared  to  the  incident  wavelength,  one  usually  obtains  the  gross  scattering 
properties  of  the  medium.  This  is  the  so-called  Rayleigh  or  low  frequency  limit, 
and  yields  corrections  to  point  scatterers.  As  far  as  the  distribution  of  the 
scatterers  is  concerned,  we  either  have  regular  arrays  of  scatterers  or  a  random 
distribution.  In  the  first  case,  we  employ  suitable  lattice  sum  while  in  the 
latter  case,  we  use  a  configuration  averaging  procedure.  If  the  scatterers  are 
sparsely  distributed,  (i.e. ,  the  concentration  is  small)  we  may  employ  a  single 
scattering  or  first  Bora  approximation. 

Most  of  the  previous  computational  results  using  scattering  theory  are,  however, 
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often  confined  to  the  Rayleigh  or  low  frequency  limit,  a  specific  angle  of 
incidence,  single  scattering  approximation  and  simple  geometries  such  as  circular 
cylinders  and  spheres.  Our  methodology  which  employs  the  T-matrix  approach  and 
suitable  statistical  averaging  procedure,  is  completely  general  and  can  be  used 
for  multiple  scattering  analyses  of  acoustic,  electromagnetic  and  elastic  waves 
emphasizing  the  inherent  unity  of  the  approach.  The  analysis  provides  a  powerful 
computational  scheme  for  determining  the  coherent  field  characteristics  of  a 
medium  containing  arbitrary' shaped  scatterers  over  a  wide  range  of  frequencies . 

We  consider  N  number  of  discrete  scatterers  of  arbitrary  shape  either  randomly 
or  periodically  distributed.  Both  parallel  and  random  orientations  of  the 
scatterers  are  considered.  When  N  is  very  large  (such  that  the  number  of  scatte¬ 
rers  per  unit  volume  is  finite),  the  T-matrix  formalism  is  combined  with  a  statis¬ 
tical  analysis  to  provide  a  self-consistent  multiple  scattering  theory.  For  a 
random  distribution  of  scatterers  we  use  a  configurational  averaging  procedure, 
while  for  a  periodic  distribution  we  use, a  suitable  lattice  sum  based  on  crysta¬ 
llographic  theory.  Inhomogeneous  media  in  which  scatterers  have  a  specific 
orientation  exhibit  anisotropy  and  polarization  effects.  These  effects  are, 
however,  nullified  for  the  case  of  randomly  oriented  scatterers.  The  information 
about  the  orientation  of  the  scatterers  has  been  incorporated  (in  terms  of  a 
rotation  matrix  consisting  of  Euler  angles)  into  the  T-matrix  of  the  scatterer 
itself  thus  making  formalism  a  convenient  computational  scheme  to  study  the 
anisotropic  effects  in  an  inhomogeneous  medium.  The  statistically  averaged 
equations  are  then  solved  using  Lax's  quasicrystalline  approximation  to  yield 
the  propagation  characteristics  of  the  coherent  waves  in  the  medium.  Analytical 
results  are  obtained  for  the  dispersion  relation  in  the  Rayleigh  or  low-frequency 
limit  for  both  2-  and  3-dimensional  scatterer  geometries.  Numerical  results  are 
presented  for  propagation  speeds,  attenuation  and  frequency  dependent  properties 
at  higher  frequencies  to  demonstrate  the  broad  applicability  of  the  T-matrix 
approach.  The  approach  presented  in  this  article  has  been  successfully  employed 
to  many  practical  scattering  problems  on  acoustic,  electromagnetic  and  elastic 
wave  fields  (22-29) . 


II.  T-MATRIX  FORMULATION  OF  MULTIPLE  SCATTERING 


We  consider  N  number  of  arbitrary  shaped  scatterers  with  a  smooth  surface  S  which 
are  referred  to  a  coordinate  system  as  shown  in  Fig.  1.  The  scatterers  may  either 
be  randomly  distributed  or  periodically  arranged.  The  orientation  of  the  scatte¬ 
rers  may  be  quite  general.  Here,  we  consider  both  specific  and  random  orientations 
of  the  scatterers.  In  Fig.  1,  0.  and  0.  refer  to  the  centers  of  the  i-th  and  j-th 
scatterers,  respectively  and  they  are  J referred  to  the  origin  0  by  the  spherical 
polar  coordinates  (rt,  8j, ♦«)•  .For  two  dimensional  problems,  we  use  cylindrical 
polar  coordinates  (rj,  ej )  .  r^  refers  to  the  vector  connecting  0i  and  0..  P  is 

any  point  in  the  medium  outside  the  scatterers  (which  we  call  the  matrix  medium) . 

We  now  describe  the  medium  and  the  scatterers  for  all  three  wave  fields,  For an 
acoustic  problem,  we  consider  fluid  scatterers  immersed  in  mother  fluid,  bubbles 
in  a  fluid,  elastic  or  viscoelastic  scatterers  immersed  in  a  fluid,  etc.  For  an 
electromagnetic  scattering  problem,  we  consider  dielectric  scatterers  in  "ee 
space,  dielectric  scatterers  embedded  in  a  different  dielectric  medium,  etc.  For 
an*alastic  wave  scattering  probl...  we  consider  elastic  or  viscoelastic  inclusions 
embedded  in  another  elastic  or  viscoelastic  material,  stress  free  or  fluid  filled 
cavities  and  cracks  in  an  elastic  or  viscoelastic  material,  etc.  The  properties  of 
the  medium  and  the  scatterers  are  given  in  terms  of  Lame* constants  X , u  and  density 
p  for  an  elastic  material,  compressibility  Xf  and  density  pf  for  non-viscous  fluids 
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Fig.  1.  Randomly  distributed  and  Fig.  2.  Translation  of  the  coordinate 

oriented  scatterers  system  of  the  j-ch  and  1-th 

scatterers . 


and  relative  dielectric  constant  e  and  permeability  u  with  respect  to  free  space 
describing  dielectric  medium,  fc'e  use  subscript  1  to  r  denote  these  qualities 
inside  the  scatterers. 

A  time  harmonic  plane  wave  of  unit  amplitude  and  frequency  to  is  incident  on  the 
medium  such  that  the  direction  of  propagation  of  the  incident  waves  is  along  the 
z-axis.  The  incident  wave  field  may  then  be  represented  by 


iCkpz 


(2.1) 


where  k^  and  are  the  congressional  (longitudinal)  and  transverse  (shear)  wave 

numbers,  respectively,  and  t  is  the  time.  The  acoustic  waves  are  purely  compre- 
ssional  type  and  thus,  the  second  term  of  Eq.  (2.1)  is  set  equal  to  zero;  for 
electromagnetic  waves  which  are  transverse  type,  the  first  term  on  the  right  hand 
side  of  (2.1)  is  zero;  for  elastic  waves  which  contain  both  compressional  and 
transverse. types,  all  the  term  of  (2.1)  are  present.  For  acoustic  and  elastic  wave 
problems,  51°  refers  to  the  incident  displacement  field  vector,  while  for  electro¬ 
magnetic  case,  it  refers  to  the  incident  electric  field  vector.  In  Eq.  (2.1),  we 
use  the  superscript  (o)  to  indicate  an  incident  wave.  The  wave  numbers  k_  and  k$ 
are  given  by 


w/c  ;  k  *  oi/c  , 
p  s  s 


(2.2) 


respectively,  where 


c  ■  /  (X  ♦  2u)/p 

P 


(elastic  waves) 


(2.3) 


/ o , 


(acoustic  waves  in  fluids) 


(2.4) 
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Cs  *  /u  /  0 

(elastic  waves) 

(2.5) 

-  C//  v  e 

(electromagnetic  waves) 

(2.6) 

In  Eqs.  (2.3)  -  (2.6),  Cp  and  cs  refer  to  compress ional  and  transverse  wave  velo¬ 
cities,  respectively.  In  Eq.(6),  C  refers  to  the  velocity  of  light  in  free  space. 
The  corresponding  quantities  inside  the  scatterers  are  differentiated  by  subscript 
1.  For  brevity,  we  use  the  notation  kT  for  the  wave  numbers  ;  r  ■  1  corresponds  to 
congressional  wave  while  r  *  2,3  corresponds  to  shear  wave. 

When  the  wave  impinges  on  a  scatterer,  part  of  the  incident  wave  is  scattered  back 
into  the  matrix  (medium  outside  the  scatterer)  and  the  rest  is  refracted  into  the 
scatterers  and  transmitted  back  into  the  matrix. We  represent  the  scattered  field  by 

u* (r)  and  the  total  field  inside  the  scatterers  by  ut (r) .  Since  the  incident  wave 
has  time  dependence  given  by  exp(-iwt),  all  these  field  quantities  have  the  same 
time  factor  which  will  be  suppressed  henceforth.  Both  these  fields  and  the. inci¬ 
dent  wave  field  satisfy  the  vector  Helmholtz  equation,  see  for  example  (30) .  Even 
though  for  acoustic  wave  propagation  in  fluids  one  could  work  with  velocity  poten¬ 
tial  and  scalar  Helmholtz  equation,  we  prefer  to  use  the  vector  displacement  vector 
field  and  the  vector  Helmholtz  equation  for  the  reasons  mentioned  in  our  previous 
papers  (20,31,32).  The  problem  at  hand  reduces  to  computing  the  total  wave  field 
at  any  point  in  the  matrix  satisfying  the  appropriate  boundary  condition  on  the 
surface  of  the  scatterers  and  radiation  conditions  at  infinity. 

The  total  field  at  any  point  in  the  matrix  (outside  the  scatterers)  is  the  sum  of 
the  incident  field  and  the  fields  scattered  by  all  the  scatterers.  This  is  written 
as 


u(r)  ■  u° (r)  ♦  £  u’(r  -  r  )  (2.7) 

i-1  1  1 

where  u*(r  -  rj)  is  the  field  scattered  by  the  i-th  scatterer  to  the  point  of  obser¬ 
vation  r.  The  field  that  excites  the  i-th  scatterer  is  the  incident  field  u°  plus 
the  fields  scattered  from  all  other  scatterers.  The  exciting- field  term  3*  is  used 
to  distinguish  between  the  field  actually  incident  on  a  scatterer  and^the  external 
incident  field  u°  produced  by  a  source  at  infinity.  Thus,  at  a  point  r  in  the 
vicinity  of  the  i-th  scatterer,  we  write 

u?(r)  *  u°(r)  u!(r  -  r.)  ,  a  <  |r  -  r  |  s  2a  (2.8) 

1  ji»i  1  J  1 

where  ’a'  is  the  radius  of  the  imaginary  sphere  circumscribing  a  scatterer.  In 
this  analysis,  we  have  assumed  that  there  is  no  interpenetration  of  the  imaginary 
spheres  of  radius  'a'  which  circumscribe  each  scatterer. 

The  T-matrix  formalism  of  scattering  we  employ  here  is  based  on  the  extended 
boundary  condition  method  which  has  been  discussed  by  many  authors  in  this  book 
(33) .  This  formalism  differs  from  the  eigen  function  expansion  technique  in  that 
the  same  basis  sets  namely,  the  vector  spherical  wave  functions  for  3-D  problems 
and  the  vector  cylindrical  wave  functions  for  2-D  problems,  may  be  used  for 
scatterers  of  any  closed  boundary  S  with  a  continuous  turning  normal.  The  vector 
fields  are  expanded  in  terms  of  a  complete  set  of  basis  functions  which  form 
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solutions  to  the  vector  Helmholtz  equation  and  are  given  by 


where 


$ 

loam 

(r)  - 

♦ 

(r)  * 

lomn 

A 

♦ 

(r)  • 

5  ■ 

[. 

mn 

ln 

r  C: 

nmn  “ 

with 

en  * 

Jh(kr,P"(co  s6,(C05m‘i”I'\l  (2.9) 

^  P  y  “  [  r.  P  "  \,ini«;c>  2 /  J 

[.  /cos  1  \1 

r  h  (k  r)  E,  (cos*)  (  I  (2.10) 

n  s  *  \sin  m*;a*  2  / J 


C2.ll) 


t(2n  *1)  (n  -  m)  ! 
n  4ir(n  +  m) ! 


(2.12) 


For  brevity,  we  abbreviate  these  vector  basis  functions  as 

*lcmn  ;  *2amn  1  *3omn  *  1  T  *  1'2’3*  For  acoustic  wave  Problems,  we  have 

only  $ln;  for  electromagnetic  waves,  we  have  <?2n  and  ^3n  l  for  elastic  waves,  we 
have  all  three  of  them.  In  Eqs.  (2.9  -  2.11),  we  have  used  spherical  polar  co¬ 
ordinates  r,0,$  with  the  origin  of  the  coordinate  system  inside  S,  hfl(l)  are  the 
spherical  Hankel  functions  of  order  n,  p“  are  the  associated  Legendre  polynomials 
and  m  is  an  integer  that  takes  values  0,  1,  2,  . . . ,  n  ;  n  *  0,  1,  2,  ...,«  for  a •  1 
and  n  *  1,2,3,  fora*  2.  Field  quantities  that  are  regular  at  the  origin 

are  expanded  in  terms  of  the  regular  (Re)  basis  set  (Re  w tan)  obtained  by  replacing 
h_  in  the  above  equations  by  jn,  the  spherical  Bessel  functions  of  the  first  kind 
of  order  n.  For  two  dimensional  problems,  we  employ  cylindrical  basis  functions, 
see,  for  example  (31) . 

We  now  expand  the  exciting  and  scattered  fields  in  terms  of  basis  functions  Re 

and  ?  ,  respectively  : 

tn 


tkCr) 


•  n 

x  n*0  Z* 


2  i  A 

? ,  Re  ^ 


(r  -  r.) 


b  _  Re  _ .  s 

tn  ml 


i  ->i 
bm  Re  4tn 


•*S 

u^.  (r) 


B*Cj)  V  BJ  V 

xn l  xn  ^  rn  Tn  n 


(2.13) 

(2.14) 


where  the  superscripts  i  or  j  on  the  basis.functions  refer  to  expansions  with 
respect  to  0,  and  0,,  respectively,  and  B-.  and  b‘n  are  unknown  coefficients  to 
be  evaluated.  The  J choice  of  basis  set  in  (2.14)  satisfies  the  radiation  condi¬ 
tion  at  infinity  for  the  scattered  field,  while  the  choice  in  (2.13)  satisfies  the 
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regularity  of  the  exciting  field  in  the  region  0< 
Substituting  (2.13)  and  (2.14)  in  (2.8),  we  obtain 


It  can  be  seen  that  the  second  series  on  the  right-hand  side  of  (2. IS)  are  expre¬ 
ssed  with  respect  to  the  center  of  the  j-th  scatterer.  In  order  to  express  them 
with  respect  to  the  i-th  center,  we  employ  the  addition  theorems  of  vector  spheri 
cal  harmonics  for  3-D  problems  and  of  Bessel  functions  for  2-D  problems.  This 
translation  has  been  discussed  in  detail  in  (34,35)  and  employed  for  acoustic, 
elctromagnetic  and  elastic  wave  problems  (22-29).  For  the  sake  of  uniformity,  we 
reproduce  the  essential  equations  which  translate  the  basis  function  from  j-th 
center  to  the  i-th  center  : 


Expressions  for  aC£,v,q)»  b(£,v»q),  a(n,£|-u,  v|q)  and  a(n,  ^-u,v|q,q-l)  are  given 
by  Cruzan  (34)  and  Fig.  2  depicts  the  geometry  of  the  translation  between  the  j-th 
and  i-th  scatterers.  In  Eq.  (2.16),  for  acoustic  wa^e  scattering,  we  keep  only 
;  for  electromagnetic  wave  scattering,  we  keep  $2  and  i-j  ;  for  elastic  wave 
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scattering,  we  need  all  three  of  then.  For  two  dimensional  problems,  the  transla¬ 
tion  is  given  by 

«  *  (r-rj)  -  C-l)n  Zd-Dm  Re  *  (r-ri)  ♦.  C^-rj)  (2.18) 
m  to  J  „  tm  T,n-m  J 

which  contain  cylindrical  Hankel  and  Bessel  functions.  For  details,  we  refer  to 
(22,  26,  36).  For  the  sake  of  simplicity  we  use  an  abbreviation  oTn  T-n-  to 
represent  the  translation  of  the  basis  functions  from  the  j-th  center  to  the  i-th 
center  and  then  we  write 


•♦■j  -*■  -r  -t-i 

r  •  (x-rj  «  2-  Cri”rj)  R« 

m  tn  J  x'n'  tn  t  n  A  ■>  T  n 


(2.20) 


It  then  remains  to  expand  the  incident  wave  also  in  the  form  of  a  series  centered 
at  the  i-th  scatterer.  Referring  to  Fig.  1,  we  can  rewrite  Eq.(2.1)  with  exp(-iwt) 
term  suppressed 


-0  .  .  .“."I*1!) 


u  ■  e 


(2.20) 


Expanding  the  terms  exp(ik„zi)  and  exp(iksZi)  in  (2.20)  in  terms  of  spherical 
Bessel  functions  and  Legendre  polynomials,  and  using  the  integral  representation 
for  spherical  Bessel  function  of  the  first  kind  and  orthogonality  relations,  we 
can  express  (2.20)  in  the  following  form,  see  for  example  Stratton  (37)  : 


i  kp 


_L_  £  £  «*  *  «  i  *•  *lt!  ‘.,0 

i  k  s-0  t-  -s  AT:s 


*TT  • 


lkp'i  ■  i 


y  T  2S+-1  i3  |  Re  A  fs  +s(s+l)6  ; 

s-1  t»-s  S(S+1)  |  21:5  L  r'l 


-r  *‘4,  [\.i 


(2. 21) 


where  6_  is  the  Kronecker  5.  For  two  dimensional  problems,  the  incident  waves 
are  expanded  by  using  Fourier  series  expansion  in  complex  form,  see  for  example 
(22,26).  For  the  sake  of  simplicity,  we  write  the  incident  wave  field  in  terms 
of  expansion  co-efficients  arn  as  follows  i—i  .  . 


-*-i 
Re  t? 

tn 


i  k  .r. 

T  1 


&  a  yU 

g|  ssl- 

jo  ft  o  Vl  i 
“J  tn  q  -h  / 

S  U8s4  i  6 


(2^  22) 


J  "  M 


Multiple  Scattering  of  Acoustic,  Electromagnetic  5  Elastic  Waves 


In  writing  (2.22),  we  have  used  the  concept  that  for  k  -  kz,  kz*r  ■  kz*(r  -  r^  ♦ 
kz*^  and  z^r.*?,.  For  acoustic  wave  scattering  (t  ■  1) ,  a^  is  equal  to  the 
first  term  of  Eq.  (2.21)  multiplying  exp  (i  k^  and  Re  ;  for  electromag¬ 
netic  waves  (t  ■  2,3),  aTp  is  equal  to  the  second  term  of  Eq.  (2.21)  multiplying 

exp  (ike)  and  Re  iT  •  Since  the  incident  wave  field  is  a  known  field,  a  are, 
r  si  Tts 

hence,  known  expansion  field  coefficients.  Substituting  Eqs.  (2.19)  and  (2.21)  in 
(2.15)  and  using  orthogonality  of  basis  functions,  we  obtain  the  following  relation 
between  the  unknown  expansion  coefficients  bTn  of  the  exciting  field  and  the 
unknown  expansion  coefficients  B  of  the  scattered  field  : 


It  has  been  shown  by  the  previous  papers  on  T-matrix  approach  that  if  the  total 
field  outside  a  scatterer  is  the  sum  of  the  incident  and  the  scattered  fields, 
the  unknown  scattererd  field  expansion  coefficients  can  be  related  to  the  incident 
field  expansion  coefficients  through  the  transition  or  T-matrix.  We  extend  this 
definition  to  our  present  problem.  Since  (tt|  +  )  is  the  total  field  at  any 

point  in  the  matrix,  the  expansion  coefficients  of  the  field  scattererd  by  the 
i-th  scatterer  may  be  formally  related  to  the  coefficients  of  the  field  exciting 
the  i-th  scatterer  through  the  T-matrix  : 


In  its  expanded  version,  Eq.  (2.24)  can  be  written  as 


The  elements  of  the  T-matrix,  T  involve  surface  integrals  which  can  be 

rn,r  n 

evaluated  in  closed  form  for  cylindrical  geometry  (2-D)  and  spherical  geometry 
(3-D),  while  for  obstacles  of  arbitrary  shape  they  can  be  evaluated  numerically 
The  T-matrix  for  a  single  scatterer  is  of  the  form 


where  Q(Re,  Re)  and  Q(0u,  Re)  are  matrices  which  are  functions  of  the  surface  S 
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of  the  scatterer  and  of  the  nature  of  the  boundary  conditions .  This  form  is  quite 
common  for  acoustic,  electromagnetic  and  elastic  wave  problems,  except  when  there 
are  solid-fluid  interfaces.  A  detailed  derivation  of  the  T-matrix  for  a  solid  in 
a  fluid  and  fluid  in  a  solid  can  be  found  in  the  paper  by  Varadan  (31) . 

With  the  scattered  field  coefficients  B^  expressed  in  terms  of  exciting  field 
coefficients  bl  and  the  T-matrix  as  given  by  (2.24),  Eq.  (2.23)  gives  the  exci¬ 
ting  field  formulation  of  the  multiple  scattering.  If  we  multiply  both  sides  of 
Eq.  (2.23)  by  the  T-matrix,  then  we  obtain  the  scattered  field  formulation  of 
multiple  scattering-  which  may  be  written  as 

B™  5  B™°  •  £  Tin,,"n"  [Vn"  ¥  > 


£  £  Bt'n'  t"  n"  (ri'rj'  ] 

1H  T'n' 


(2.27) 


In  its  expanded  version,  Eq.  (2.27)  can  be  written  as 

Ls™l  fcT11,^  (Ti2,s(i)  ai3,r,: 
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I 

I 

I 
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P 


In  Eq.  (2.29),  denotes  the  sum  over  all  scatterers  except  the  i-th.  Thus,  we 
have  eliminated  the  unknown  exciting  field  expansion  coefficients  through  the  use 
of  the  T-matrix  resulting  in  a  set  of  equations  involving  the  expansion  coeffi¬ 
cients  of  the  scattered  field  only.  If  the  scatterers  are  all  identical,  then 

*  T.  If  they  are  of  different  sizes,  we  introduce  a  suitable  size  distri¬ 
bution  for  the  scatterers  and  find  average  value  of  the  T-matrix. 

It  can  be  seen  from  Eq.  (2.27)  that  the  scattererd  field  coefficients  explicitly 
depend  on  the  position  and  orientation  of  the  scatterers.  Depending  on  what  infor¬ 
mation  we  choose  to  put  into  the  distribution  function,  we  can  make  our  model 
more  realistic  but  the  analysis  in  turn  gets  more  complicated.  Here,  we  consider 
both  random  and  periodic  distributions  with  specific  and  random  orientations. 


III.  RANDOM  DISTRIBUTION  OF  SCATTERERS  WITH  SPECIFIC  ORIENTATION 


We  consider  N  number  of  randomly  distributed  homogeneous  scatterers.  The  scatte¬ 
rers  are  assumed  to  have  a  specific  orientation  (with  major  axis  parallel  to  the 
y-axis)  as  shown  in  Fig.  3. 


Distribution  functions  are  given  by  the  total  probability  density  of  finding  the 
scatterers  at  positions  r^,  Tj,  ...»  r^  etc.,  which  in  turn  can  be  written 
in  terms  of  conditional  probabilities  : 


PC^)  p(rx  r2. 


p(*ji*i)  pCV 


(3.1) 
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In  Eq.  (3.1),  p(r.)  denote  the  probability  of  finding  a  scatterer  r.  and  p(r.|  r.) 

*  ^  1  J  1 

denotes  the  conditional  probability  of  locating  a  scatterer  at  r.  given  that  a 
scatterer  is  located  at  r  etc.  A  prime  in  the  first  of  Eq.  (24)  means  is 
absent  while  two  primes  in  the  second  of  Eq.  (24)  means  both  r.  and  r.  are 
absent.  We  also  denote  the  configurational  averaging  of  a  statistical  ^quantity 
f  as 


<f>i  -/•••  f  f  PC?,.  r2, 

•'V  ■'V 

<f>ij  afm“  J  f  p<^i»  '2* 


• • • ’ »  |  r .)  d?1  d?2  . . . ' ,  ....  d?N 


•••’’  •”  ?N  '  V  V  d^ld*2  ’  *  *  *  * 


...  drN  (3.2) 


Equation  (3.2a)  implies  that  we  have  averaged  over  all  scatterers  except  the  i-th, 
while  Eq.  (3.2b)  implies  that  we  average  over  all  scatterers  except  the  i-th  and 
j-th,  and  so  on. 


If  the  scatterers  are  randomly  distributed,  the  positions  of  all  scatterers  are 
equally  probable  within  volume  V,  and  hence  its  distribution  is  uniform  with 
probability  density 


V 


V 


(3.3) 


where  nQ  is  the  uniform  number  density  of  the  scatterers  and  N  is  the  total 
number  of  scatterers.  In  addition,  for  non-overlap  of  the  imaginary  spheres  cir¬ 
cumscribing  each  scatterer,  we  approximate  the  conditional  probability  as  follows  : 


l^i  *  rj|  >  2a 
l*i  -  rj|  <  2a 


(3.4) 


where  'a'  is  the  radius  of  the  circumscribing  sphere.  A  suitable  correlation  in 
the  position  may  also  be  added  to  (3.4)  but  is  omitted  here  for  simplicity.  The 
form  of  pair  correlation  in  Eq.  (3.4)  describes  the  usual  radially  symmetric 
distribution  function  with  an  exclusion  surface  or  "hole"  corresponding  to  a 
sphere  of  radius  2a.  Other  types  of  exclusion  surfaces  are  discussed  in  recent 
papers  by  Twersky  (8,  9).  We  refer  the  reader  to  Twersky's  comments  on  choosing 
radially  symmetric  statistics  for  non-symmetrical  scatterers,  namely  making  the 
region  of  non-overlap  of  two  scatterers  to  be  a  circle  or  sphere  circumscribing 
the  scatterers. 

Multiplying  both  sides  of  Eq.  (2.2,6)  by  the  probability  density  given  by  Eq.  (3.1) 
and  using  (3.2)  -  (3.4),  we  obtain  the  configurational  average  of  the  scattered 
field  coefficients  : 

i  ^  r  i  k  .r 

<Btn  ”i  *  L,..  Tm,T"n"  a-r”n"  e 
T"n'  L 
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♦  Ju- 
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E  E  [<>;- n->i: 

t'n'y 
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aT'n',T"n"  j 


(3.5) 


where  V'  denotes  the  volume  of  the  medium  excluding  the  volume  of  a  sphere  of 

radius  2a.  For  identical  scatterers,^  can  be  replaced  by  (N-l) .  The  above 

equation  indicates  that  the  configurational  average  with  one  scatterer  fixed  is 
given  in  terms  of  the  configurational  average  with  two  scatterers  fixed.  We  choose 
to  solve  these  equations  by  using  Lax's  quasicrystalline  approximation  where  we 
assume  that  there  is  no  correlation  between  the  scatterers  except  that  of  non¬ 
overlapping  which  can  be  written  mathematically  as 


<  B  >. .  *  <  B  _>i 
rn  i]  tn  j 


(3.6) 


Equation  (3.6)  also  implies  that  the  neighbourhood  of  every  scatterer  is  the 
same.  For  identical  scatterers,  Eq.  (3.5)  may  be  rewritten  using  Eq.  (3.6)  as 
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In  its  expanded  version,  Eq.  (3.7)  can  be  written  as 
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(3.9) 


This  is  a  system  of  integral  equations  for  the  unknown  coefficients  <  B*n  >i- 
Similar  expressions  for  the  average  exciting  field  coefficients  may  be 
obtained  from  the  exciting  field  formalism  of  multiple  scattering. 


116 


V.K.  Varadan 


IV.  PROPAGATION  CHARACTERISTICS  OF  THE  AVERAGE  WAVES  IN  THE  MEDIUM 


To  solve  the  integral  equations  given  by  (3.7),  we  consider  the  inhomogeneous 
medium  with  discrete  scatterers  as  a  homogeneous  continuum  and  assume  that  the 
average  coherent  wave  is  a  plane  wave  propagating  with  an  effective  wave  number 
K  in  the  same  direction  as  the  incident  plane  wave.  We  can  thus  write 


„i  ' 

<  B  >  »  X 

TO  TO 


(4.1) 


where  XTn  is  the  amplitude  of  the  coherent  wave. 

Substituting  Eq.  (4.1)  in  (3.7),  employing  divergence  theorem  to  convert  the 
volume  integral  in  (3.7)  to  surface  integrals  and  using  the  extinction  theorem 
which  cancels  the  incident  wave,  see  for  example  (8,9,22),  we  obtain  a  set  of 
simultaneous  coupled  homogeneous  equations  for  the  coefficients  X^  given  by 

ln'+n"l  q  JHo 

5|n--ni  V"'  {*.-  **)«* 

(4.2) 

where  c  *  4ir  a3  nn/3  is  the  effective  spherical  concentration  of  the  scatterers 
per  unit  volume,  0  C^  is  an  expression  containing  Wigner  coefficients,  and 


JHq  -  2  kT.a  j (2  k'a)h'(2  kT>a)  -  2  K  a  hq(2  kT-a)  jq  (2  Ka)  (4.3) 

Equation  (4.2)  can  be  written  in  its  expanded  form  as  follows  : 
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a(ml'n1l”P'ml^ 


CT12^ 


b(n1,m,q) 


a(m^,nj|-p,m|q, 


(4.4a) 

(4.4b) 

(4.4c) 


Eauation  (4.4b)  can  be  obtained  from  (4.4a)  by  replacing  T1*,  T^.T^byT^1,  T^, 
T‘3,  while  (4.4c)  can  be  obtained  by  replacing  T11,  T12,  T13  by  T31,  T32,  and  T33. 

a 

For  acoustic  wave  problem,  we  get  uncoupled  equation  for  X  '  ;  for  electromagnetic 

wave  problem,  we  get  coupled  equations  in  terms  of  X*  and  n  X*  ;  for  elastic  wave 

o  zn  a  3n  a 

problem,  we  obtain  coupled  equations  in  terms  of  Xjn,  X2n  and  X|p.  Similar 
equations  can  also  be  obtained  from  the  average  exciting  field  coefficients,  if 
one  chooses  to  use  exciting  field  formalism. 


Equation  (4.2)  is  a  system  of  simultaneous  linear  homogeneous  equations  for  the 
unknown  amplitudes  X  .  For  a  nontrivial  solution,  we  require  that  the  determinant 
of  the  truncated  coefficient  matrix  vanishes,  which  yields  an  equation  for  the 
effective  wave  number  K  in  terms  of  kT  and  the  T-matrix  of  the  scatterer.  This  is 
the  dispersion  relation  for  the  scatterer  filled  medium.  Equation  (32)  is  a  gene¬ 
ral  expression  valied  for  any  arbitrary  shaped  scatterer,  since  the  T-matrix  is 
the  only  factor  that  contains  information  about  the  exact  shape  and  boundary 
conditions  at  the  scatterer.  Thus  the  formalism  presented  here  is  valid  for  all 
the  three  wave  fields.  The  effective  wave  number  K  obtained  in  the  analysis  is  a 
complex  quantity,  the  real  part  of  which  relates  to  the  phase  velocity,  while 
the  imaginary  part  relates  to  attenuation  of  coherent  waves  in  the  medium. 


V.  PERIODIC  DISTRIBUTION  WITH  PARALLEL  ORIENTATION 


For  periodic  distribution  the  analysis  will  introduce  the  geometry  of  the  lattice 
leading  to  different  results  for  different  packing,  and  is  similar  to  the  one 
employed  in  crystallographic  theories,  see  for  example  (38,39).  Here  we  apply 
lattice  sum  and  include  contributions  from  nearest  neighbors.  One  couls  use  other 
refinements  following  the  theories  presented  in  (8,9,38,39).  The  averaging  over 
the  position  of  individual  scatterers  can  be  easily  accomplished  since  there  is 
no  restriction  on  the  position  of  one  scatterer. 

A  procedure  for  evaluating  a  lattice  sum  over  a  simple  lattice  is  to  obtain  by 
direct  summation  the  contributions  by  the  lattice  points  within  a  certain  radius  Rj 
given  by 


3  V 


and  to  replace  the  summation  over  the  points  beyond  this  radiusjjy  an  integral.  In 
Eq.  (S.l),  Nq  is  the  number  of  scatterers  counting  the  point  rj  "  0  and  V  is  the 
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volume  of  the  unit  cell,  (see  (40)  for  a  given  typ«  of  an  array)'.  The  analysis  is 
presented  in  (41, 42)  and  the  resulting  scattered  field  coefficients  are  given  by 


Z  T1 

T"n"  Tn,T"n” 


>  V  ri 
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t  n  ,r"n' 


(5.2) 


To  obtain  a  solution  to  the  above  equation,  we  assume  a  plane  wave  propagating 
with  an  effective  wave  number  K  in  the  same  direction  as  the  incident  wave  : 


tn 


x  .ilc-r‘ 

tn 


(5.3) 


Substituting  Eq.  (S.3)  in  Eq.  (S.2)  and  following  the  procedure  we  had  outlined 
in  the  previous  section,  we  obtain  a  set  of  simultaneous  coupled  homogeneous 
equations  for  the  coefficients  XTfl  given  by 
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(5.5) 


and  c  is  the  concentration  as  defined  before. 


VI.  ARBITRARY  ANGLE  OF  INCIDENCE  AND  RANDOM  ORIENTATION 


Inhomogeneous  media  in  which  scatterers  have  a  specific  orientation  will  exhibit 
anisotropy  and  polarization  effects  when  the  waves  are  incident  at  arbitrary  angles 
to  the  symmetry  axes  of  the  scatterers.  When  the  incident  wave  direction  is  along 
the  symmetry  axis  of  the  oblate  spheroidal  scatterers,  we  found  that  the  effect 
on  polarization  is  zero,  and  the  coherent  wave  propagates  with  an  effective  wave 
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number  K  as  if  it  propagates  in  an  effective  homogeneous  and  isotropic  continuum, 
see  for  example  (25) .  When  the  scatterers  are  randomly  oriented,  the  effects  of 
anisotropy  and  polarization  are  also  nullified  (29) . 

For  synnnetrically  oriented  scatterers,  we  define  the  Euler's  angles  a,S,y  of  the 
symmetry  ares  of  the  scatterers  (x,y,z)  with  respect  to  the  laboratory  co-ordinate 
system  (X,Y,Z),  see  Fig. 4.  All  quantities  that  are  referred  to  the  x,y,z  system 
are  distinguished .by  a  a.  Then,  one  could  write  the  spherical  harmonics  using  the 
relation  operator -as  follows  (for  details,  we  refer  to  (43))  : 


Y 


Lao 


L 


Dmm'  oo' 


(o,3,y) 


(6.1) 


where  D  is  the  rotation  matrix  associated  with  rotation  operator.  The  rotation 
operator  leaves  the  length  of  the  position  vector  |r|  invariant.  The  rotation 
matrix  can  be  easily  incorporated  into  the  T-matrix  to  obtain  a  new  T-matrix  which 
can  be  written  as 


T  «  (DV1  T  (D*) 


(6.2) 


where  D*  is  the  matrix  transpose  of  D.  Equation  (6.2)  gives  the  desired  relation 
between  tfc*  T-matrices  evaluated  with  respect  to  the  two  set  of  coordinate  axes. 

T  is  independent  of  position  and  orientation  and  is  hence  the  same  for  identical 
scatterers.  The  matrices  T,  however,  is  different  if  the  orientation  of  the 
scatterers  is  not  the  same.  Substituting  Eq.  (6.2)  in  Eq.  (4.2)  and  Eq.  (5.4),  we 
obtain  the  dispersion  relation,  phase  velocity,  coherent  wave  attenuation  and 
polarization  effects  for  an  arbitrary  angle  of  incidence  for  random  and  periodic 
distribution,  respectively.  This  is  one  of  the  basic  advantage  of  formulating  the 
multiple  scattering  in  terms  of  a  T-matrix. 


For  random  orientation,  one  has  to  average  over  all  Euler  angles  a,0,y,  and  the 
information  can  thus  be  incorporated  into  the  T-matrix  : 


<T>  ■ 


D(a,6,r) 


T  D‘1(a,8,Y) 


sin  6  d6  dy  da  • 


(6.3) 


With  Eq.  (6.3),  Eqs.  (4.2)  and  (5.4)  yield  the  dispersion  relation  for  both 
random  and  periodic  distributions, respectively,  with  random  orientation. 


VII.  RESULTS  AND  CONCLUSIONS 


Using  the  theories  outlines  in  previous  sections,  we  present  some  analytical  and 
numerical  results  for  variety  of  2-D  and  3-D  problems  in  all  three  wave  fields 
to  show  the  broad  applicability  of  our  multiple  scattering  approach.  The  sample 
examples  given  may  find  applications  in  many  fields  of  engineering  and  science 
as  outlined  in  the  introduction.  We  present  results  for  two  dimensional  parallel 
cylinders  of  elliptical  cross  section  randomly  and  periodically  (hexagonal  array) 
distributed  with  specific  and  random  orientation.  The  wave  is  incident  normal  to 
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the  cylinders.  The  polarization  of  the  wave  is  along  the  axis  of  the  cylinder 
(SH-waves).  We  also  consider  spherical,  oblate  spheroidal  scatterers  of  various 
aspect  ratios  randomly  distributed  with  specific  and  random  orientations.  At  low 
frequencies,  analytical  expressions  are  derived  for  the  effective  wave  number  in 
the  average  medium  as  a  function  of  the  geometry ,  the  material  properties  and  the 
angle  of  orientation  of  the  scatterers.  The  formulation  is  ideally  suited  for 
numerical  computation  of  phase  velocity  and  attenuation  at  higher  frequencies  as 
evidenced  by  the  results  presented  here  and  in  (22  -  29) .  In  addition,  we  present 
numerical  results. for  dynamic  shear  moduli  for  2-D  case  as  a  function  of  frequency 
using  the  T-matrix  approach  and  the  work  of  Bedeux  and  Mazur  (44)  and  Varadan  and 
Vezzetti  (45) . 


Vila.  Rayleigh  or  Low  Frequency  Limit 

In  the  Rayleigh  or  low  frequency  limit,  the  size  of  the  scatterers  is  considered 
to  be  small  when  compared  to  the  incident  wavelength.  It  is  then  sufficient  to 
take  only  the  lowest  order  coefficient  in  the  expansion  of  the  fields.  In  this 
limit,  the  elements  of  the  T-matrix  can  be  obtained  in  closed  form  for  various 
simple  shapes  (46).  It  can  be  shown  that  at  low  frequencies,  only  XtQ,  XTland  X^. 

of  Eq.  (4.2)  or  Eq.  (5.4)  make  a  contribution.  After  some  manipulations  of  the 
resulting  3*3  determinant,  we  obtain  the  dispersion  relations  : 

(i)  Parallel  Elliptic  Cylinders  (Elastic  SH-Waves) 

a)  Specific  Orientation  of  Cross  Section  (see  Fig.  S  and  paper  (23)) 


Fig.  5.  Random  distribution  of  aligned  scatterers  -  oblique  incidence 
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and  c 


where 


Cii) 


■  (2ic/ir2)  a^^  e 


■i2a 


1  -  2k2a2  tnfk^ajJ 


e2 


ir  -  4ic  Sq  8ic  SQ  k#2  a2  tn(k3a) 


(1/ir2)  ir  -  4ic  SQ  *  8ic  SQ  k52  a2  in(k3a)J 
x  ir  -  2ic  |  1  -  2ks2a2  in(ksa) 


(4c/ir)  a^  sin  2  o  £l  -  2ks2  a2  In(k^a) 


CSX  ♦  &1  e‘i2a)] 


•  Cl/ir)  [*  -  2ic  {l  -  2ks2a2  ta(ksa)|  (Sx  -  ax  e'i2o)j 

■  (2ic/ir2)  Sx  [it  ♦  8ic  SQ  ks2a2  tn(ksa)] 

■  (1/ir2)  £ir  ♦  8ic  Sq  k*  a2  in(k3a)J  .  £ir  2ic  CS^-*-  aj  ei2°)J 

•  i(ir/8)  (d  -  1)  (b/a)  ^2  -  (d  -  1)  ks2ab  ln(ksa)]  (7.2) 


il  S'— *.*&<*>  f ,  +  n  .  (i)  (1  .  ffl)  rCi+b2/a2) 

8  (b/a  ♦  m)  (1  *  mb/a)  l  4  L 

x  (  1  ♦  m2)  4  mb/aj  ks2ab  Zn(k3a)/(1  +  mb/a)  (b/a  +  m)| 

fl.-  alfl  -  faj/alUh/al  fi  .  m  1  fS  ♦  b/a,)2Xl.^a 
(b/a  ♦  m) (1  ♦  mb/a)  L  4  (b/a  *  m) (1  +  mb, 

ks2ab  in(ksa)J 


JJL 

8 


-  ”2lv 
mb/ a) 


ira2nQ  is  the  concentration  of  the  circumscribing  circle. 


b)  Random  Orientation 


K2/k2  -  (1  ♦  Cjtj)  [l  ♦  ct(d  -  1)]  /(I  -  cjtj) 


(7.3) 


tj  ■  (1  -  m2) (a  ♦  b)2  /  4(b  ♦  ma)  (mb  ♦  a) 


(7.4) 


Random  Distribution  of  2-D  cracks  (Elastic  SH  -  Waves) 
a)  Parallel  Orientation 


(7.S) 


K2  (4  »  ira2nQ)2-  (ira2n0)2  e~**°  -2i(ira2n0)2  e  *2asin  2<x 


ks2  jd  -  ira2nQ(l  -  e12o)J  ^4  ♦  ira2nQ(l  ♦  e*i2a)|-  2i(ira2nQ)25in  2a 
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b)  Random  Orientation 


K2/ks2  ■  (1  *  ira2nQ/4)  (1  -  ira2nQ/4) 


(7.6) 


where  2a  is  the  length  of  the  crack  and  nQ  is  its  number  density,  and  a  is  the 
angle  the  incident  wave  makes  with  the  major  axis.  In  these  equations,  d  -ratio 
of  densities  of  the  scatterer  to  that  of  the  matrix  -  pj/p  and  m  ■  ratio  of  shear 
moduli  of  the  scatterer  to  that  of  the  matrix  -  y./y  and  a  and  b  are  semi  major 
and  minor  axes  of  the  elliptic  cylinder,  respectively.  The  analysis  for  both 
compressional  and  shear  waves  incident  normal  to  cylinders  is  presented  in  (27) . 


(iii)  Spheres  (elastic  waves) 


•(f)' 


(1  ♦  9c  Et)  (1  ♦  3c  Eq)  [l  ♦  3c  -fe-  (2  +  ^)] 


1  -  ISc  E2  1  +  3c  Eq  ♦  c  E2  (2  ♦  ) 

Lb  • 


Hi1" 

V 


(7.7) 


where 


(f)‘ 


(1  <•  9C  E,H1  *  -f  «  E,  ^  ] 


1  ♦  -3- 


t  ‘  V4  ‘  ^ 


(7.8) 


E0  “ 


E1  " 


1  3X  +  2y  -  (3Xi  ♦  2y-i) 
3  4y  ♦  $Xj*  2y^  U” 


4-  (4-  - « 

9  P 


y  (y^y)  24u1(u1-u)  -  (Xi*2yi)  (Wy^tty) 

24y1(y1-y)  -  (X1*2y1)  (19Ul«-  16y)  • 


(7.9) 


■4y(Uj-y)  ♦  3(X  ♦  2y)  (2y^  ♦  3y) 

and  c  -  4  ira3nQ/3  is  the  concentration  of  spheres. 

(iv)  Spheres  in  Free  Space  (Electromagnetic  Waves), 

1  ♦  2c  ISl 


H)  ■ 


_r-i  *± 


crx  ♦  2 


1  -  c 


(7.10) 
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Spheroids  in  Free  Space  CElectTomagnetic  Waves 


Random  Orientation 


where  c  is  the  concentration  of  the  circumscribing  spheres  and  fj,  f2  and  f3 
are  functions  of  the  eccentricity  of  the  spheroids  e  *  [(a/b)2  -  l]1*  and  are 
given  in  (25) . 


In  the  Rayleigh  limit,  the  value  of  K  as  determined  by  the  above  dispersion  rela¬ 
tions  is  a  real  quantity  for  lossless  material  and  a  complex  quantity  for  lossy 
material,  and  relates  to  phase  velocity  V_  ■  w/K.  In  this  limit,  we  normally  study 
the  dependence  of  phase  velocity  on  concentration,  angle  of  incidence  and  aspect 
ratio  of  the  scatterers.  In  Figs.  6  and  7,  we  have  plotted  the  normalized  phase 
velocity  Vn/cs  versus  concentration  for  various  values  of  angle  of  incidence  a 
for  2-D  (SH-Wave)  cylinders  randomly  distributed  with  parallel  and  random  orien¬ 
tations.  We  assume  the  density  ratio  d  »  2.53/2.72  and  shear  modulus  ratio 
m  ■  2S/3.87  (no  loss  in  the  material)  for  obtaining  results  in  Fig.  6  while  Fig.  7 
gives  the  result  for  cracks.  The  general  tendency  of  the  phase  velocity  is  to 
increase  (for  inclusions)  and  decrease  (for  cracks  and  cavities)  as  concentration 
increases.  The  results  also  indicate  that  phase  velocity  decreases  as  the  angle 
of  incidence  a  increases.  For  o»0  and  a»ir/2,  our  results  agree  with  well  known 
results  for  cracks  (3).  It  can  easily  be  seen  from  Eq.  (7.2)  that  for  a  given 
concentration  and  a,  the  phase  velocity  decreases  with  decreasing  b/a.  In  Figs. 

6  and  7,  we  have  also  plotted  the  corresponding  results  for  random  orientation. 

Equation  (7.10)  is  recognized  as  the  dispersion  relation  of  the  Clausius  -  Mcssotti 
form.  The  dispersion  relation  for  spheroidal  scatterers  given  by  Eq.  (7.11) 
appears  to  be  new  ;  it  reduces  to  Eq.  (7.10)  in  the  limit  e  •+•().  In  Fig.  8,  we  have 
plotted  the  normalized  phase  velocity  V_/cs  for  dielectric  spheres  and  spheroidal 
scatterers  in  free  space  randomly  distributed  and  oriented  versus  concentration. 
Both  real  and  complex  dielectric  constants  (e_),  of  the  scatterers  are  assumed. 


The  dispersion  relation  given  in  Eqs.  (7.1)  and  (7.3)  may  also  be  useful  in 
obtaining  the  effective  shear  modulus  at  low  frequencies.  Defining  - 


uotaming  me  exxeccive  snear  modulus  at  low  trequencies.  Defining  an  average 
shear  modulus  <y>  ■  w 2  <p>  JC2  where  <p>  «  c  (1  -  c)p  is  the  average  density 
we  find  the  effective  shear  modulus  of  the  composite  medium  containing  cylindri¬ 
cal,  inclusions  : 
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Fig.  8a.  Normalized  phase  velocity  vs.  concentration  for  spherical  and  spheroidal 
dielectric  scatterers  in  free  space;  -  lossy,  -  lossless 


Fig.  8b.  Attenuation  vs.  concentration  c  for  spheroidal  lossy  dielectric 
scatterers  in  free  space 


<u>/u  •  (<p>/p)  (fcs/K)2 
for  more  detail,  we  refer  to  (28). 

Using  the  dispersion  relation  given  in  Eqs.  (7.7)  and  (7.8),  one  could  also 


(7.12) 
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compute  the  effective  Shear  and  Bulk  moduli  of  an  elastic  material  containing  a 
random  distribution  of  stress  free  bubbles  or  cavities  : 


<U>  _  4u  -  3c  Et (9A  ♦  14u)  .  <B>  C5X  ♦  2y)  [1  -  6c  Enl  (7.13) 

U  4u  +  6c  E2(3X  >  8u)  '  B  (3A  *  2u)  [1+  3c  Egj 

where  <u>  and  <B>  are  the  effective  Shear  and  Bulk  moduli  of  the  inhomogeneous 
medium,  and  u  and  B  are  the  Shear  and  Bulk  moduli  of  the  matrix  medium.  This 
result  might  find  -applications  in  geophysics  and  material  science,  see  for  example 
Chaban  (48) . 


Vllb.  Phase  Velocity  and  Attenuation  at  Higher  Frequencies 


To  study  the  response  at  resonant  and  higher  frequencies,  we  must  consider  higher 
powers  of  kTa,  and  this  implies  that  a  larger  number  of  terms  (XT_)  nfust  be  kept 
in  the  expansion  of  the  average  field.  This  is  best  done  numerically.  For  a  given 
value  of  ka,  the  T-matrix  for  the  scatterer  is  computed.  Next,  the  coefficient 
matrix  M  corresponding  to  (Eqs.  (4.2)  and  (5.4))  is  formed.  The  complex 

determinant  of  the  coefficient  matrix  is  computed  using  standard  Gauss  elimination 
techniques.  For  a  given  kTa,  the  root  of  the  equation  det  M  =  0  is  searched  in 
the  complex  K  plane  (K^  *  iK^)  using  Muller’s  method.  Good  initial  guesses  were 
provided  by  the  Rayleigh  limit  solutions  at  low  values  of  kTa  and  these  could  be 
used  systematically  to  obtain  convergence  of  roots  at  increasingly  higher  values 
of  kTa.  The  real  part  Kj  determines  the  phase  velocity,  while  the  imaginary  part 
K2  determines  the  coherent  wave  attenuation.  Here,  we  present  results  for  2  -  D 
SH-Wave  and  3-D  electromagnetic  problems.  We  define  the  normalized  phase  velo¬ 
city  as  Vp/cs  *  ks/Kj  and  the  attenuation  coefficient  ■  4ir  K2/K’  so  as  to  make 
it  dimensionless. 

For  two  dimensional  elliptical  cylinders  (SH-Waves) ,  phase  velocity  and  attenua¬ 
tion  are  plotted  versus  ksa  in  Figs.  9-13  for  various  values  of  b/a,  o  and  actual 
concentration  c^  *  nab  ng  for  both  random  distribution  with  specific  and  random 
orientations  and  periodic  distribution  with  specific  orientation.  We  assume 
d  *  2.53/2.72  and  m  =  25/3.87.  The  general  behavior  of  the  phase  velocity  is  to 
increase  at  relatively  low  frequencies  and  then  decrease  rather  sharply  for  higher 
frequencies.  For  b/a  =  1.0,  the  results  agree  qualitatively  with  those  obtained 
by  Sutherland  and  Lingle  (49).  As  a  increases  from  0?  to  90®,  it  has  been  found 
that  the  normalized  phase  velocity  is  reduced  as  evidenced  by  Fig.  9.  The  aniso¬ 
tropic  effect  is  greatly  pronounced  for  lower  values  of  b/a  indicating  that 
there  is  a  substantial  reduction  in  phase  velocity  for  nearly  flat  oriented 
inclusions  when  the  incident  wave  is  in  the  direction  perpendicular  to  the 
inclusions.  This  observation  is  of  importance  in  geophysics.  The  results  in  Fig.  9 
indicate  that  there  is  a  significant  difference  in  phase  velocity  for  various 
orientations  and  that  the  corresponding  values  for  random  orientation  lie  in 
between  0®  and  90®  orientations. 

In  Figs.  10  and  12,  we  have  plotted  the  coefficient  of  attenuation  S<j  versus  ksa. 
These  results  indicate  that  anisotropic  effect  is  greatly  pronounced  ;  substantial 
reduction  of  attenuation  is  observed  when  o  »  90®.  The  general  behavior  of  the 
damping  or  attenuation  is  to  increase  rapidly  with  frequency  at  low  frequency 
range  and  decrease  sharply  or  fluctuates  slightly  when  frequency  is  increased. 

The  results  also  indicate  that  there  is  a  tendency  for  the  attenuation  to  increase 
rapidly  again  at  higher  frequencies.  The  attenuation  calculations  for  random 
orientation  lie  between  0®  and  90®  orientations,  see  Fig.  10.  In  Fig.  10,  we 
have  also  plotted  the  coefficient  of  attenuation  versus  ksa  for  2  -  D  cracks 
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I.  Aligned  scatterers 
a  “  90° 

II.  Aligned  scatterers 

V  o-60° 

\  III.  Aligned  scatterers 
\  o  ■  0° 

\  IV.  Randomly  oriented 
^  \  -  scatterers 


^  C^-0.5 


/  IV 


'V — £ _ -KH 
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II.  Aligned  scatterers^ 
/  <x-90° 

/  II.  Aligned  scatterers 
o-60° 

III.  Aligned  scatterers 
o-0° 

IV.  Randomly  oriented 
scatterers 

V.  Random  distribution 
of  cracks 


Fig.  9.  Normalized  phase  velocity  vs.  k  a  Fig.  10.  Attenuation  S  vs.  kaa  for  a 
for  a  random  distribution  of  3  random  distribution  of 

aligned  and  randomly  oriented  aligned  and  randomly  orient- 

scatterers  ed  scatterers  and  cracks . 

i  60 1 - ; - ; - ’ - 1 


Fig.  11.  Normalized  phase 
velocity  vs.  k  a 
for  a-45°  and 
C,«0.5 


Multiple  Scattering  of  Acoustic,  Electromagnetic  $  Elastic  Waves 


129 


randomly  distributed. 


The  frequency  at  which  Sd  begins  to  decrease  is  usually  referred  to  as  cut  off 
frequency  of  the  first  pass  band.  At  the  cut  off  frequency,  Sd  decreases  because 
energy  begins  to  pass  into  the  second  pass  band,  see  the  explanations  in  (49,50). 
For  periodic  distribution,  these  cut  off  frequencies  seem  to  move  towards  lower 
ksa  as  the  intercylindrical  spacing  between  the  cylinders  increase,  see  Fig.  13. 


In  Figs.  14  -  17,  .-we  have  plotted  phase  velocity  and  attenuation  coefficients  for 
spherical  and  spheroidal  dielectric  scatterers  in  free  space  randomly  distributed 
with  parallel  and  random  orientations.  The  calculations  were  performed  for  both 
lossless  and  lossy  dielectric  scatterers.  For  lossless  scatterers,  we  assume 
dielectric  constant  er  *  3.168  which  corresponds  to  ice  at  microwave  frequencies. 
The  imaginary  part  of  the  dielectric  constant  for  ice  is  relatively  small  when 
compared  to  the  real  part,  see  Ray  (47).  For  lossy  case,  we  consider  rain  parti¬ 
cles  which  have  complex  dielectric  constants  given  as  functions  of  temperature 
and  frequency.  The  complex  dielectric  constants  of  rain  particles  are  taken  from 
the  paper  by  Ray  (47) .  For  our  calculations ,  we  assume  the  temperature  to  be 
equal  to  5*C. 


Fig.  14a.  Normalized  phase  velocity 
vs  k  a  for  aligned 
spheroidal  dielectric 

scatterers 


Fig.  14b.  Normalized  phase  velocity 

vs  k  a  for  randomly  oriented 
spheroidal  dielectric 
scatterers 
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0  0*5  1*0  1*5  2 


♦Fig.  17 

Fig.  15.  Attenuation  Sd  vs  lcga  for  lossy 

dielectric  scatterers  for  various 
values  of  a/b  and  for  c»10“^ 

Fig.  16.  Attenuation  Sd  vs  k,a  for  lossless 
aligned  dielectric  scatterers  for 
various  values  of  c  and  a/b*1.25 

Fig.  17  Attenuation  vs  k#a  for  lossless 
dielectric  scatterers  for  various 
values  of  a/b  and  for  c*0.2 
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In  Fig.  IS,  we  show  the  attenuation  coefficient  for  lossy  dielectric  scatterers 
as  a  function  of  ksa  for  concentration  c  ■  10~4.  The  results  are  plotted  for 
spherical  and  oblate  spheroidal  scatterers  randomly  distributed  with  both  parallel 
and  random  orientations.  For  the  case  of  parallel  orientation,  the  spheroids 
are  assumed  to  be  oriented  with  their  axes  of  revolution  along  the  z-axis  (see 
Fig.  3),  and  the  incident  wave  is  also  assumed  to  propagate  along  the  z-axis.  In 
both  specific  and  random  orientation  cases,  it  is  well  known  that  the  bulk  medium 
is  not  anisotropic  and  is  characterized  by  a  single  wave  number  K.  However,  when 
the  wave  is  incident  obliquely  to  parallelly  oriented  scatterers,  we  obtain  the 
anisotropic  effects  (just  like  the  2  -  D  case  discussed  above)  and  the  bulk  medium 
will  then  be  characterized  by  an  effective  wave  number  with  two  different  polari¬ 
zations  . 

In  Figs.  16  and  17,  we  have  plotted  Sd  for  lossless  spheroidal  dielectric  scatte¬ 
rers  for  two  different  concentrations,  c  »  0.1  and  0.2,  and  for  two  different 
aspect  ratios,  a/b  »  1.2S,  2.  The  significant  feature  in  the  case  of  c  ■  0.2  is 
the  presence  of  a  sharp  null  at  ksa  »  0.75. 

The  multiple  scattering  analysis  presented  so  far  can  also  be  extended  to  compute 
effective  dynamic  moduli  of  a  composite  medium  (28) .  Bedeaux  and  Mazur  (44)  have 
obtained  expressions  for  the  average  dielectric  tensor  in  the  medium  and 
Varatharajuiu  (Varadan)  and  Vezzetti  (45)  have  shown  that  for  a  given  statistical 
model  the  zeros  of  the  propagator  in  this  model  medium  yield  the  dielectric 
tensor  defined  by  Bedeaux  and  Mazur  (44) .  In  (28) ,  we  have  extended  this  discu¬ 
ssion  to  elastic  wave  propagation  wherein  we  have  given  sample  calculations  for 
dynamic  shear  moduli  as  a  function  of  ksa. 

In  conclusion,  we  have  presented  a  multiple  scattering  formalism  for  coherent 
wave  propagation  of  acoustic,  electromagnetic  and  elastic  waves  through  an 
inhomogeneous  medium  composed  of  either  random  or  periodic  distribution  of 
2-dimensional  and  3-dimensional  scatterers.  The  scatterers  are  assumed  to  have 
either  specific  or  random  orientation.  An  important  advantage  is  realized  through 
the  use  of  the  T-matrix  to  characterize  the  scattering  properties  of  any  one 
scatterer.  The  methodology  adapted  in  this  analysis  is  general  and  can  be  used  to 
include  such  effects  as  scatterer  size  distribution,  depolarization,  etc* 
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